Phase space tweezers for tailoring cavity fields by quantum Zeno dynamics 
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We discuss an implementation of Quantum Zeno Dynamics in a Cavity Quantum Electrodynamics 
experiment. By performing repeated unitary operations on atoms coupled to the field, we restrict 
the field evolution in chosen subspaces of the total Hilbert space. This procedure leads to promising 
methods for tailoring non-classical states. We propose to realize 'tweezers' picking a coherent field at 
a point in phase space and moving it towards an arbitrary final position without affecting other non- 
overlapping coherent components. These effects could be observed with a state-of-the-art apparatus. 

PACS numbers: 03.65.Xp, 42.50.Dv, 42.50.Pq 
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In the Quantum Zeno effect (QZE) [1], repeated pro- 
jective measurements block the evolution of a system in 
a non-degenerate eigenstate of the measured observable. 
It has been observed on two-level systems [2] and on an 
harmonic oscillator in a Cavity Quantum Electrodynam- 
ics (CQED) experiment [3|. A Quantum Zeno dynamics 
(QZD) [4] takes place when the system is not confined 
to a single state, but rather evolves under the action of 
its free hamiltonian i7 in a multidimensional subspace 
of its Hilbert space. This can be achieved either by re- 
peated measurements of an observable with degenerate 
eigenvalues, or by repeated actions of a unitary kick Uk 
with multi-dimensional invariant subspaces, the two pro- 
cedures being physically equivalent [5]. We focus here 
on the latter case, related to the so-called 'bang-bang' 
control [G"] and NMR manipulation techniques [7]. 

The system evolution is stroboscopic, alternating small 
free evolution steps described by U{5t) = ex.-p{—iHSt/h) 
with Uk kicks. The succession of N steps (fixed du- 
ration t = Ndt) corresponds to the unitary Uz{N) = 
[UKU{t/N)]^. It is, in the TV ^ oo limit, the evolution 
under an effective Hamiltonian Hz = XI u P^iHP^ where 
the ^'s are the projectors on the invariant subspaces of 
Uk [J]- By choosing properly Uk (or equivalently the 
repeatedly measured observable), one can tailor the sys- 
tem evolution, leading to decoherence control [5|, state 
purification [8] and quantum gates implementation [9|. 
QZD can also inhibit entanglement between subsystems, 
making a quantum evolution semi-classical [10]. 

In this Letter, we propose a CQED implementation of 
the bang-bang QZD control. We exploit the non-linearity 
of the atom-cavity system [ll[ to implement a photon- 
number-selective Uk- The field dynamics in its Hilbert 
space 1-L is confined in two disconnected subsets H<s and 
H>s, corresponding to photon numbers smaller or larger 
than a preset value s. This leads to novel methods of 
non-classical field states preparation and tailoring. We 
propose a 'phase space tweezer' picking selectively a co- 




FIG. 1. (a) Proposed experimental scheme. A slow atomic 
beam extracted from a 2D-Magneto Optical Trap (bottom) 
forms an atomic fountain with atoms nearly at rest in the 
center of the high-quality microwave Fabry- Perot cavity C 
(only one mirror shown). Sources S and S' address respec- 
tively the dressed atomic levels and the cavity mode. Elec- 
trodes around the cavity mirrors generate the electric fields 
preparing the circular state shown in the center, (b) Scheme 
of the dressed atomic levels. The arrow indicates the photon- 
number selective transition addressed by S for s = 1. 



herent state component of a quantum superposition and 
moving it at will in phase space independently from the 
other components. 

Our proposal could be implemented in a microwave 
CQED experiment with circular Rydberg atoms and a su- 
perconducting millimeter- wave cavity [11| . The cavity C 
[Fig. HJa)] is crossed by a slow beam of Rubidium ground 
state atoms, in an 'atomic fountain' arrangement. Close 
to their turning point, atoms are nearly at rest. One of 
them is promoted to the circular level h (principal quan- 
tum number 49) using static and r.f. electric fields. This 
operation does not change the state of C, tuned on reso- 
nance with the 51.1 GHz transition between g and e (50 
and 51 circular states). The analysis that follows pertains 
to the dynamics of this single atom and the cavity. 

The source S drives the h ^ g transition near 
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FIG. 2. (a) QZD dynamics in l-i<Q. Ten snapshots of the field Wigner function W{^) obtained after a number of steps indicated 
above each frame. The cavity is initially in its vacuum state, s = 6 and /3 = 0.1. The EC is plotted as a blue dashed line, (b) 
QZD dynamics in 1-L>q. Same as (a) with an initial a = — 5 amplitude, (c) Same as (b), with an initial amplitude a — — 4 + i\/6- 
In (b) and (c) the successive frames correspond to the same step numbers as in (a). 



54.3 GHz. It does not couple directly to the off-resonant 
cavity. It realizes Ur by probing the eigenstates of the 
atom/cavity system [Fig. [IJb)]. The energies of the |/i, n) 
states (atom in h with n photons) do not depend on the 
atom-cavity coupling, since C is far off-resonance from 
the h ^ g transition. The source 5* couples \h^n) to the 
dressed states |±,n) = (|e,n — 1) ± |^,n))/>/2, which 
are superpositions of the degenerate uncoupled states 
|e,n — 1), \g^n). The splitting between dressed states 
is Qy/n where Q is the vacuum Rabi frequency. Hence, 
the \h,n) -> |+,n) transition frequency depends upon 
n. We tune S to perform on the \h^s) -^ |+,5) tran- 
sition a 27r Rabi pulse whose amplitude is weak enough 
(and its duration correspondingly long enough) not to 
appreciably affect \h^n) with n ^ s. It results in the 
transformation |/i,n) -^ (— l)'^'^^|/i,n). The atom always 
ends up in h while the field experiences \Jk = Us with 
Us = t — 2\s){s\. Such a photon- number dependent Rabi 
pulse [12| was used with 5 = 1 for a single-photon QND 
detection ^ and for a CNOT gate in CQED [14]. 

The free cavity dynamics is produced by the source 
S\ resonantly coupled with C and acting during time 
intervals 5t between two Us operations. Being not reso- 
nant with the atom in /i, S^ leaves it unaffected. The 
free evolution is described by the Hamiltonian H = 
—i{£*a — £a^)^ where £ is the source amplitude and a 
{a^) the photon annihilation (creation) operator. We use 
an interaction representation eliminating the field phase 
rotation at cavity frequency. The unitary U{5t) is the dis- 
placement D{p) = exp(/3at - /3*a), with (3 = SSt/h <C 1. 
After p repetitions of UsU{dt)^ the cavity state can be 
reconstructed [15]. 

The eigenvalues of Us are —1 and +1. The former cor- 
responds to the one-dimension eigenspace Hs, generated 
by \s) (projector P^). The latter is associated to the di- 



rect sum of H<s (projector P<s), generated by the Fock 
states |0), . . . , |5 — 1), and Hys (projector Pys) generated 
by Fock states above \s). The projectors P^ (/i = +, — ) 
are thus P_ = P^ and P+ = P<s + P>s- 

Since H is a. linear combination of a and a^ ^ Hz reduces 
to P^sHP^s + PysHPys- Under the QZD, field states 
restricted to H<s and H>s remain confined in these sub- 
spaces, \s) realizing a hard 'wall' between them. This wall 
induces remarkable features in the evolution. If we start 
from the vacuum in C with 5 = 1, the system remains 
inside H<i i.e. in |0). We recover the QZEj3i]- 

We have simulated this QZD procedure [16|. Fig. [Jfa) 
presents 10 snapshots of the field Wigner function, W{^)^ 
separated by intervals of 5 steps, for 5 = 6 and j3 = 0.1. 
The field starts from |0) G H<6. Its amplitude first in- 
creases along the real axis (free dynamics). Between 15 
and 20 steps, the amplitude reaches :^ 2 and QZD comes 
into play. The coherent state 'collides' on the t/g-induced 
'wall', materialized in phase space as an 'Exclusion Cir- 
cle' (EC) of radius ^6 (dashed line in Fig. EJ. The field 
amplitude stops growing and undergoes a progressive tt 
phase shift between steps 20 and 30. At step 25, the field 
is in a 'cat state', quantum superposition of two compo- 
nents with opposite phases. The fringing feature inside 
the EC is the signature of the quantum coherence of this 
superposition. This cat contains only odd photon num- 
bers (probabilities for 5, 3 and 1 photons are 0.63, 0.31 
and 0.03). At step 35, the field state is nearly coherent 
with an amplitude close to -2. It then resumes its mo- 
tion from left to right along the real axis, going through 
|0) again (around step 45) and heading towards its next 
'collision' with the EC. 

Fig. [3] presents the long-term evolution of the field 
energy. During the first few hundred steps, the oscil- 
lations reveal the quasi-periodic motion of the field in- 
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FIG. 3. Energy of the field as a function of tiie number of 
steps for a QZD dynamics inside 7i<c6- Conditions are the 
same as for Fig. [2ja), which presents snapshots of the first 
oscillation. 



side the EC. State distortions, however, accumulate and 
damp the oscillations, whose contrast nearly vanishes af- 
ter 800 steps. Since there is only a finite set of frequencies 
in P<6i^P<6, we observe at longer times a quantum re- 
vival [ll|. The energy oscillations resume and the field 
comes back to an oscillating coherent state periodically 
colliding with the EC as described above. 

Fig. Efb) illustrates QZD in H>s with snapshots of the 
field Wigner function for s = 6 and an initial coherent 
state \a = —5). The field collides on the EC after 20 
steps. It undergoes a QZD-induced tt phase shift being, 
after 25 steps, in a cat state. After 30 steps, the state is 
again nearly coherent with a positive amplitude. It re- 
sumes its motion along the real axis. After 45 steps, its 
amplitude is slightly larger than 4.5. It would be —0.5 in 
the case of free dynamics. The QZD-induced phase inver- 
sion accelerates the 'propagation' in phase space. This 
opens interesting possibilities when the initial amplitude 
is such that the field state collides tangentially on the 
EC. The parts of the Wigner function that come clos- 
est to the EC propagate faster than others. The state is 
distorted and ends up strongly squeezed [Fig. Efc)]. 

QZD can be generalized to ECs centered at an arbi- 
trary point 7 in phase space by changing the kick op- 
erator Uk from Us to /7s (7) = D{^)UsD{—^) (these 
displacements 7 being also performed by S'). After 
p steps, the global evolution operator is Uz{s^j^p) = 
[Us{j)D{/3)]P which can be expressed, using displace- 
ment operator commutation relations, as Uz{s^^^p) — 
D{j)Uz{s, 0,p)D{-j) exp[2ip^(/37*)]. Up to a topologi- 
cal phase, the state after p steps is equivalently obtained 
by first displacing the field by —7, performing p QZD 
steps in an EC centered at origin and finally displacing 
back the field by 7. 

This leads to the concept of phase space tweezer. Ap- 
plying this procedure with 5 = 1 to an initial cat state 
1 7) + I Of) {{jla) ^ 0), we can selectively block the evo- 
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FIG. 4. (a)-(b) Initial and final Wigner functions W(^) for 
a phase space tweezer operation. The first step ECs are de- 
picted as solid lines in (a) and dotted lines in (b), the final 
ECs as solid lines in (b). The arrows in (b) indicate the two 
EC centers trajectories, (c) Cat created by a vacuum state 
crush operation. Initial (dotted lines) and final (solid line) 
ECs are plotted, with arrows indicating the motion of their 
centers, (d) Four-component cat created by three successive 
crushing operations. 



lution of 1 7) while leaving the other component free 
to evolve. After N steps, we get the 'stretched' cat 
\j)^D{Nf3)\a) = |7)-fexp[zAr^(/3a*)]|a + Ar/3). 

In an interesting variant, the position 7^ of the tweezer 
at step p is changed by a small amount at each step 
(|7p+i — 7p| <C 1), while /3 is set to zero (no free evo- 
lution). The sequence of {7^} defines a trajectory T in 
phase space followed by the center of the EC. A coherent 
state with the initial amplitude a^ = 71 follows adiabat- 
ically this trajectory and becomes, after the p-th step, a 
coherent state with amplitude ap = 7p. We realize in this 
way a tweezer which moves one selected coherent state, 
while not affecting the evolution of all the coherent states 
whose amplitude remains away from T. 

A combination of tweezers can move all the com- 
ponents of a superposition of non-overlapping coherent 
states from arbitrary initial to final positions. An obvi- 
ous method is to grasp them one by one, driving them 
from their initial to their final position (taking care to 
move them so that different components never overlap). 
The tweezers can also be used in parallel by applying 
incremental motions alternatively on each component. 

Fig. |4] illustrates this procedure for a two-component 
cat, initially \a) -\- \ — a) with a = 2. It is turned in 
100 steps (50 on each component) into \a') + | — a') with 
a' = 5i. Panels (a) and (b) present the initial and final 
Wigner functions. The final fidelity is 98.8% with respect 
to the expected cat. It remains greater than 68% if the 



operation is performed in 20 steps only, exhibiting the ro- 
bustness of this adiabatic procedure, which is promising 
for experimental implementations. 

The initial cat state can be generated in various ways 
using dispersive [17| or resonant [18] atoms. An adapted 
tweezer procedure also prepares from |0) mesoscopic su- 
perpositions which are approximations of these cat states. 
We 'crush' the vacuum between two s = 1 ECs whose 
centers, initially at ±2.5, move simultaneously towards 
each other in 200 steps, until they reach the origin (the 
wavefunction remaining outside both ECs). Fig. HJ^c) 
presents the final Wigner function. The state is a super- 
position of two well-separated components (average en- 
ergy 6.4 photons). Its fidelity with respect to a |a) + |— a) 
cat with the same average energy is 42%. The final com- 
ponents can be crushed again and so on, leading to a 
superposition state with an arbitrary number of com- 
ponents. Fig. m^c) presents the Wigner function of a 
four-component cat obtained by crushing again each of 
the two components in Fig. lU^b) between ECs moving 
towards each other along the imaginary axis direction. 

QZD is also obtained when the interrogation pulse has 
a Rabi angle different from 27r. The pulse performs 
then a unitary kick acting on the atom/cavity system, 
which mixes \h,s) with |+,5) and would create atom- 
field entanglement if C contained s photons. This unitary 
admits an invariant subspace, belonging to the eigenvalue 
-hi spanned by the projection \h){h\ {P<s + P>s)^ the 
same as for a 27r pulse. Starting from an atom in \h) and 
a field in H<s or H>s, we obtain a QZD leaving the atom 
in \h) and the field in its initial subspace. Under perfect 
QZD, the cavity never contains s photons and the atom 
and field are never entangled by the interrogation pulse. 
Obviously, QZD is not achieved if is very small, each 
kick operation being too close to 1. We have checked 
numerically that, for ^ c^ 1, we recover within a good 
approximation all the results described above. 

Numerical simulations can also take into account re- 
alistic experimental imperfections. We have simulated a 
typical microwave CQED experiment [15] with a cavity 
damping time Tc = 0.13 s and an atom-field coupling 
Q/27T = 50 kHz. We have taken into account the limited 
selectivity of the interrogation pulse, which must be at 
the same time much shorter than Tc and much longer 
than 1/Q. We have carefully optimized the interroga- 
tion pulse parameters [19] for a tweezer operation lead- 
ing, in 10 steps, from an |a = —2) -\- \a = 2) cat to 
\a = —3) -\- \a = 3) . The total duration of the simulated 
experiment is 3.4 ms and the final fidelity with respect 
to the ideal cat is 76%. 

Quantum Zeno dynamics applied to a field oscillator 
leads to novel methods for tailoring non-classical fields. 
Phase space tweezers are promising tools for generat- 
ing and manipulating arbitrary superpositions of coher- 
ent states (i.e. arbitrary states superpositions) [19] and 



to study their evolution under the effect of decoher- 
ence ^] . Experimental demonstrations of these effects 
are within reach of a state-of-the-art microwave CQED 
set-up. They could also be implemented in circuit QED, 
a field in which high QTc quality factors are also real- 
ized [20|. Manipulating at will the state of a quantum 
oscillator in its phase space provides a new insight into 
the physics of mesoscopic quantum superpositions and 
the exploration of the quantum to classical boundary. 
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and DECLIC projects) and by the ANR (QUSCO- 
INCA). 



[1] B. Misra and E. C. G. Sudarshan, J. Math. Phys. 18, 756 

(1977); H. Nakazato et al, Int. J. Mod. Phys. B 10, 247 

(1996); Phys. Lett. A 217, 203 (1996); D. Home and M. 

A. B. Whitaker, Ann. Phys. (N.Y.) 258, 237 (1997); K. 

Koshino and A. Shimizu, Phys. Rep. 412, 191 (2005). 

[2] W. M. Itano, D. J. Heinzen, J. J. Bollinger and D. J. 

Wineland, Phys. Rev. A 41, 2295 (1990); B. Nagels, L. 

J. F. Hermans and P. L. Chapovsky, Phys. Rev. Lett. 79, 

3097 (1997); P. G. Kwiat et a/., Phys. Rev. Lett. 83, 4725 

(1999); C. Balzer, R. Huesmann, W. Neuhauser and P. 

Toschek, Opt. Comm. 180, 115 (2000); E. W. Streed et 

a/., Phys. Rev. Lett. 97, 260402 (2006); O. Hosten et a/.. 

Nature 439, 949 (2006). 

[3] J. Bernu et al, Phys. Rev. Lett. 101, 180402 (2008). 

[4] P. Facchi and S. Pascazio, Phys. Rev. Lett. 89, 080401 

(2002); J. Phys. A 41, 493001 (2008). 
[5] P. Facchi, D. A. Lidar and S. Pascazio, Phys. Rev. A 69, 

032314 (2004). 
[6] L. Viola and S. Lloyd, Phys. Rev. A 58, 2733 (1998). 
[7] W. A. Anderson and F. A. Nelson, J. Chem. Phys. 39, 
183 (1963); R. R. Ernst, J. Chem. Phys. 45, 3845 (1966); 
R. Freeman, S. P. Kempsell and M. H. Levitt, J. Magn. 
Reson. 35, 447 (1979); M. H. Levitt, R. Freeman and T. 
A. Frenkiel, J. Magn. Reson. 47, 328 (1982). 
[8] H. Nakazato, M. Unoki and K. Yuasa, Phys. Rev. A 70, 

012303 (2004). 
[9] X. Q. Shao et al, J. Opt. Soc. Am. 26, 2440 (2009). 
[10] R. Rossi, K. M. Fonseca Romero and M. C. Nemes, Phys. 

Lett. A 374, 158 (2009). 
[11] S. Haroche and J. M. Raimond, Exploring the quantum 

Oxford University Press (2006). 
[12] M. Franca Santos, E. Solano and R. L. de Matos Filho, 

Phys. Rev. Lett. 87, 093601 (2001). 
[13] G. Nogues et al. Nature, 400, 239 (1999). 
[14] A. Rauschenbeutel et ai., Phys. Rev. Lett. 83, 5166 

(1999). 
[15] S. Deleglise et al. Nature, 455, 510 (2008). 
[16] S. M. Tan, J. Opt. B 1, 424 (1999). 
[17] M. Brune et al. Rev. Lett. 77, 4887 (1996). 
[18] A. Auffeves et ai., Phys. Rev. Lett. 91, 230405 (2003). 
[19] J. M. Raimond et ai., in preparation. 
[20] B. R. Johnson et ai.. Nature Physics, 
DOL10.1038/NPHYS1710 (2010). 



